LINKING ALLOMETRIC SCALING THEORY WITH LIDAR REMOTE SENSING FOR IMPROVED BIOMASS ESTIMATION AND ECOSYSTEM CHARACTERIZATION

Accurate quantification of forest carbon stocks and fluxes is critical for the successful modeling and mitigation of climate change. This research focuses on forest carbon stock quantification, both in terms of testing emerging remote sensing approaches to forest carbon modeling, and examining allometric equations used to estimate biomass stocks in field plots. First, we test controversial theoretical predictions of forest allometry through the mapping of the allometric variability using field plots across the U.S. we find that there is considerable variability in forest allometry across space, largely driven by local environment and life history. However, in tall forests, allometries tend to converge toward theoretical predictions, suggesting that theory may be a useful constraint on allometry in certain forests. Second, we shift to an analysis of empirical allometries by developing an algorithm to extract individual crown information from forest systems and using it for biomass mapping and allometric equation testing. Third, we test whether individual tree structure bolsters biomass modeling capabilities in comparison to tradition, plot-aggregated LiDAR metrics. As part of this analysis we also test an allometric scaling-based approach to biomass mapping. We find that individual tree-level structure only improves biomass models when there is considerable spatial heterogeneity in the forest. Also, allometric scaling-based only worked in one study site, and failed in the other two sites because there was little or no relationship between basal area and maximum canopy height. Finally, we applied LiDAR datasets to an analysis of the effects of sample size on empirical allometry development. We found that small samples sizes tend to result in an under sampling of large stems, which yields a more linear fit than the true allometry. An assessment of the potential carbon implications of this problem yielded site-level biomass predictions with biases of 10-178%. We suggest that empirical allometric equations developed on small sample sizes, as applied in the U.S., yield potentially large errors in biomass and therefore require careful reassessment. In combination with our findings regarding the spatial variability of forest allometry, we believe that the limiting factor to forest carbon estimation is the use of allometric equations

The Use of Sun Elevation Angle for Stereogrammetric Boreal Forest Height in Open Canopies

Stereogrammetry applied to globally available high resolution spaceborne imagery (HRSI; less than 5 m spatial resolution) yields fine-scaled digital surface models (DSMs) of elevation. These DSMs may represent elevations that range from the ground to the vegetation canopy surface, are produced from stereoscopic image pairs (stereo pairs) that have a variety of acquisition characteristics, and have been coupled with lidar data of forest structure and ground surface elevation to examine forest height. This work explores surface elevations from HRSI DSMs derived from two types of acquisitions in open canopy forests. We (1) apply an automated mass-production stereogrammetry workflow to along-track HRSI stereo pairs, (2) identify multiple spatially coincident DSMs whose stereo pairs were acquired under different solar geometry, (3) vertically co-register these DSMs using coincident spaceborne lidar footprints (from ICESat-GLAS) as reference, and(4) examine differences in surface elevations between the reference lidar and the co-registered HRSI DSMs associated with two general types of acquisitions (DSM types) from different sun elevation angles. We find that these DSM types, distinguished by sun elevation angle at the time of stereo pair acquisition, are associated with different surface elevations estimated from automated stereogrammetry in open canopy forests. For DSM values with corresponding reference ground surface elevation from spaceborne lidar footprints in open canopy northern Siberian Larix forests with slopes less than10, our results show that HRSI DSM acquired with sun elevation angles greater than 35deg and less than 25deg (during snow-free conditions) produced characteristic and consistently distinct distributions of elevation differences from reference lidar. The former include DSMs of near-ground surfaces with root mean square errors less than 0.68 m relative to lidar. The latter, particularly those with angles less than 10deg, show distributions with larger differences from lidar that are associated with open canopy forests whose vegetation surface elevations are captured. Terrain aspect did not have a strong effect on the distribution of vegetation surfaces. Using the two DSM types together, the distribution of DSM-differenced heights in forests (6.0 m, sigma = 1.4 m) was consistent with the distribution of plot-level mean tree heights (6.5m, sigma = 1.2 m). We conclude that the variation in sun elevation angle at time of stereo pair acquisition can create illumination conditions conducive for capturing elevations of surfaces either near the ground or associated with vegetation canopy. Knowledge of HRSI acquisition solar geometry and snow cover can be used to understand and combine stereogrammetric surface elevation estimates to co-register rand difference overlapping DSMs, providing a means to map forest height at fine scales, resolving the vertical structure of groups of trees from spaceborne platforms in open canopy forests

Fusing simulated GEDI, ICESat-2 and NISAR data for regional aboveground biomass mapping

Accurate mapping of forest aboveground biomass (AGB) is critical for better understanding the role of forests in the global carbon cycle. NASA's current GEDI and ICESat-2 missions as well as the upcoming NISAR mission will collect synergistic data with different coverage and sensitivity to AGB. In this study, we present a multi-sensor data fusion approach leveraging the strength of each mission to produce wall-to-wall AGB maps that are more accurate and spatially comprehensive than what is achievable with any one sensor alone. Specifically, we calibrate a regional L-band radar AGB model using the sparse, simulated spaceborne lidar AGB estimates. We assess our data fusion framework using simulations of GEDI, ICESat-2 and NISAR data from airborne laser scanning (ALS) and UAVSAR data acquired over the temperate high AGB forest and complex terrain in Sonoma County, California, USA. For ICESat-2 and GEDI missions, we simulate two years of data coverage and AGB at footprint level are estimated using realistic AGB models. We compare the performance of our fusion framework when different combinations of the sparse simulated GEDI and ICEsat-2 AGB estimates are used to calibrate our regional L-band AGB models. In addition, we test our framework at Sonoma using (a) 1-ha square grid cells and (b) similarly sized irregularly shaped objects. We demonstrate that the estimated mean AGB across Sonoma is more accurately estimated using our fusion framework than using GEDI or ICESat-2 mission data alone, either with a regular grid or with irregular segments as mapping units. This research highlights methodological opportunities for fusing new and upcoming active remote sensing data streams toward improved AGB mapping through data fusion.</p

Biomass estimation from simulated GEDI, ICESat-2 and NISAR across environmental gradients in Sonoma County, California

Estimates of the magnitude and distribution of aboveground carbon in Earth's forests remain uncertain, yet knowledge of forest carbon content at a global scale is critical for forest management in support of climate mitigation. In light of this knowledge gap, several upcoming spaceborne missions aim to map forest aboveground biomass, and many new biomass products are expected from these datasets. As these new missions host different technologies, each with relative strengths and weaknesses for biomass retrieval, as well as different spatial resolutions, consistently comparing or combining biomass estimates from these new datasets will be challenging. This paper presents a demonstration of an inter-comparison of biomass estimates from simulations of three NASA missions (GEDI, ICESat-2 and NISAR) over Sonoma county in California, USA. We use a high resolution, locally calibrated airborne lidar map as our reference dataset, and emphasize the importance of considering uncertainties in both reference maps and spaceborne estimates when conducting biomass product validation. GEDI and ICESat-2 were simulated from airborne lidar point clouds, while UAVSAR's L-band backscatter was used as a proxy for NISAR. To estimate biomass for the lidar missions we used GEDI's footprint-level biomass algorithms, and also adapted these for application to ICESat-2. For UAVSAR, we developed a locally trained biomass model, calibrated against the ALS reference map. Each mission simulation was evaluated in comparison to the local reference map at its native product resolution (25 m, 100 m transect, and 1 ha) yielding RMSEs of 57%, 75%, and 89% for GEDI, NISAR, and ICESat-2 respectively. RMSE values increased for GEDI's power beam during simulated daytime conditions (64%), coverage beam during nighttime conditions (72%), and coverage beam daytime conditions (87%). We also test the application of GEDI's biomass modeling framework for estimation of biomass from ICESat-2, and find that ICESat-2 yields reasonable biomass estimates, particularly in relatively short, open canopies. Results suggest that while all three missions will produce datasets useful for biomass mapping, tall, dense canopies such as those found in Sonoma County present the greatest challenges for all three missions, while steep slopes also prove challenging for single-date SAR-based biomass retrievals. Our methods provide guidance for the inter-comparison and validation of spaceborne biomass estimates through the use of airborne lidar reference maps, and could be repeated with on-orbit estimates in any area with high quality field plot and ALS data. These methods allow for regional interpretations and filtering of multi-mission biomass estimates toward improved wall-to-wall biomass maps through data fusion.</p

Detecting Change in Forest Structure with Simulated GEDI Lidar Waveforms: A Case Study of the Hemlock Woolly Adelgid (HWA; Adelges tsugae) Infestation

The hemlock woolly adelgid (HWA; Adelges tsugae) is an invasive insect infestation that is spreading into the forests of the northeastern United States, driven by the warmer winter temperatures associated with climate change. The initial stages of this disturbance are difficult to detect with passive optical remote sensing, since the insect often causes its host species, eastern hemlock trees (Tsuga canadensis), to defoliate in the midstory and understory before showing impacts in the overstory. New active remote sensing technologies-such as the recently launched NASA Global Ecosystem Dynamics Investigation (GEDI) spaceborne lidar-can address this limitation by penetrating canopy gaps and recording lower canopy structural changes. This study explores new opportunities for monitoring the HWA infestation with airborne lidar scanning (ALS) and GEDI spaceborne lidar data. GEDI waveforms were simulated using airborne lidar datasets from an HWA-infested forest plot at the Harvard Forest ForestGEO site in central Massachusetts. Two airborne lidar instruments, the NASA G-LiHT and the NEON AOP, overflew the site in 2012 and 2016. GEDI waveforms were simulated from each airborne lidar dataset, and the change in waveform metrics from 2012 to 2016 was compared to field-derived hemlock mortality at the ForestGEO site. Hemlock plots were shown to be undergoing dynamic changes as a result of the HWA infestation, losing substantial plant area in the middle canopy, while still growing in the upper canopy. Changes in midstory plant area (PAI 11-12 m above ground) and overall canopy permeability (indicated by RH10) accounted for 60% of the variation in hemlock mortality in a logistic regression model. The robustness of these structure-condition relationships held even when simulated waveforms were treated as real GEDI data with added noise and sparse spatial coverage. These results show promise for future disturbance monitoring studies with ALS and GEDI lidar data

Local discrepancies in continental scale biomass maps: a case study over forested and non-forested landscapes in Maryland, USA

Continental-scale aboveground biomass maps are increasingly available, but their estimates vary widely, particularly at high resolution. A comprehensive understanding of map discrepancies is required to improve their effectiveness in carbon accounting and local decision-making. To this end, we compare four continental-scale maps with a recent high-resolution lidar-derived biomass map over Maryland, USA. We conduct detailed comparisons at pixel-, county-, and state-level. Spatial patterns of biomass are broadly consistent in all maps, but there are large differences at fine scales (RMSD 48.5–92.7 Mg ha−1). Discrepancies reduce with aggregation and the agreement among products improves at the county level. However, continental scale maps exhibit residual negative biases in mean (33.0–54.6 Mg ha−1) and total biomass (3.5–5.8 Tg) when compared to the high-resolution lidar biomass map. Three of the four continental scale maps reach near-perfect agreement at ~4 km and onward but do not converge with the high-resolution biomass map even at county scale. At the State level, these maps underestimate biomass by 30–80 Tg in forested and 40–50 Tg in non-forested areas. Local discrepancies in continental scale biomass maps are caused by factors including data inputs, modeling approaches, forest/non-forest definitions and time lags. There is a net underestimation over high biomass forests and non-forested areas that could impact carbon accounting at all levels. Local, high-resolution lidar-derived biomass maps provide a valuable bottom-up reference to improve the analysis and interpretation of large-scale maps produced in carbon monitoring systems.https://doi.org/10.1186/s13021-015-0030-

Statistical properties of hybrid estimators proposed for GEDI – NASA’s Global Ecosystem Dynamics Investigation

NASA’s Global Ecosystem Dynamics Investigation (GEDI) mission will collect waveform lidar data at a dense sample of ∼25 m footprints along ground tracks paralleling the orbit of the International Space Station (ISS). GEDI’s primary science deliverable will be a 1 km grid of estimated mean aboveground biomass density (Mg ha ^−1 ), covering the latitudes overflown by ISS (51.6 °S to 51.6 °N). One option for using the sample of waveforms contained within an individual grid cell to produce an estimate for that cell is hybrid inference, which explicitly incorporates both sampling design and model parameter covariance into estimates of variance around the population mean. We explored statistical properties of hybrid estimators applied in the context of GEDI, using simulations calibrated with lidar and field data from six diverse sites across the United States. We found hybrid estimators of mean biomass to be unbiased and the corresponding estimators of variance appeared to be asymptotically unbiased, with under-estimation of variance by approximately 20% when data from only two clusters (footprint tracks) were available. In our study areas, sampling error contributed more to overall estimates of variance than variability due to the model, and it was the design-based component of the variance that was the source of the variance estimator bias at small sample sizes. These results highlight the importance of maximizing GEDI’s sample size in making precise biomass estimates. Given a set of assumptions discussed here, hybrid inference provides a viable framework for estimating biomass at the scale of a 1 km grid cell while formally accounting for both variability due to the model and sampling error